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The measurement of the parity of two qubits is a primitive of quantum computing that allows creating 
deterministic entanglement. In the field of circuit quantum electrodynamics, a scheme to achieve parity 
measurement of two superconducting qubits has been proposed and analyzed under the usual rotating-wave 
approximation (RWA). We show that the same scheme can be carried over beyond this approximation, to the 
regime of ultrastrong coupling, with an improvement in the fidelity. 



I. INTRODUCTION 



II. PARITY MEASUREMENT 



The study of superconducting qubits coupled to a 
microwave resonator is known as circuit quantum 
electrodynamics (cQED) d I2). Experiments have 
demonstrated the coherent control Q, entanglement 
generation [4 ] and readout [ 5 , 6 ] of qubits with high fidelity, 
making circuit QED one of the most promising architectures 
to achieve fault tolerant quantum computing Q. 

One of the latest experimental achievements is the 
demonstration of ultrastrong coupling [9), where the 
coupling strength between the qubit and the resonator 
is of the order of their resonance frequencies (g/oj r > 
0.1). In this regime, physics is no longer captured by 
the usual rotating-wave approximation (RWA) and novel 
effects have been predicted QUI [TTJ. In the context of 
quantum information processing, stronger coupling suggests 
faster dynamics, which is a very welcome development, 
since decoherence is one of the limiting factors in the 
performances of circuit QED. However, the translation is not 
so straightforward: for instance, two-qubit gates must be 
re-designed lfT2l because the usual architectures would suffer 
a rapid degradation of fidelity in this new regime fl3l . 

In this paper, we explore the extension to the ultrastrong 
coupling regime of another basic procedure: parity 
measurement, which creates deterministic entanglement by 
extracting the overall parity of a two- or multi-qubit state. 
The existing proposals for parity measurement fT4ifT6l are 
based on the dispersive qubit readout: the qubits are strongly 
detuned from the cavity, the resonance frequency of the cavity 
is shifted depending on the state of the qubits, and the joint 
information on single ifTTl or multi-qubit (HI [18) properties 
is obtained by homodyne measurement on a probe field 
transmitted through the resonator. This dispersive readout 
remains possible in the ultrastrong coupling regime [19]. 
As a result, we find that the proposed schemes for parity 
measurement reach even higher fidelity when g/a> r increases. 

The paper is organized as follow. In Sec. [II] we first present 
the drive-Rabi model in the dispersive regime. Next, we 
extend the two-qubit parity measurement to the ultrastrong 
coupling domain. In Sec. [Ill] by using the example of 
entanglement generation we discuss the deviations between 
the RWA and exact solution in different regimes of the 



In the dispersive regime with multiple qubits, the oscillator 
frequency exhibits a shift depending on the collective states 
of the qubits. This feature allows the joint readout of the 
multi-qubit state parity. By restricting to the two-qubit 
scenario, the objective of a two-qubit parity measurement 
is to distinguish the states between two orthogonal parity 
subspaces: r H+ = span(|gg) , \ee)) and = span(|g£> , \eg)). 
A perfect parity measurement projects the state of the 
system onto either subspace, without gaining any single-qubit 
information in order to preserve the superpositions in the 
post-measurement states. We consider a cQED system 
consisting of two superconducting qubits coupled to a 
transmission line resonator with a driving field applied to the 
input port (Fig. [TJ. The corresponding Hamiltonian reads 
(h = I) H2QH 
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The first three terms are the usual Rabi model for atom-light 
interaction where oj a j is the frequency of qubit j, oj r is the 
frequency of the resonator and gj is the coupling strength 
between jth qubit and the resonator. The last term describes 
the feeding of the resonator by the classical field of amplitude 
e m at frequency cj m . |2D 

This drive-Rabi Hamiltonian, despite its simple form, is 
analytically non-trivial. Although the analytical solution of 
this Rabi Model has recently been found by Braak [ 22 ], the 
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parameters. Sec. IV summarizes our results. 



FIG. 1. (Color online) Schematic of the cQED setup. The microwave 
driving field is generated and transmitted through the resonator 
which is coupled ultrastrongly with the qubits. In the homodyne 
measurement, the microwave drive is amplified and mixed with a 
local oscillator (LO), where the in-phase (/) and the quadrature (Q) 
components of the transmitted field are measured. 
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general formula is in term of power series of a transcendental 
function. Since we are operating in the dispersive regime 
where the resonator and qubits frequencies are far detuned 
from each other compared to the coupling strength gj, i.e. 
gj/Aj <sc 1 and gj/^j <^ 1, where Aj = oj r - oj a j and 
X 7 = oj r + CL) a j, a good approximate solution can be obtained 
by using the transformation fT9ll23ll24l 



U = exp 
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where Xj = acr+ - a^crj, Yj - acr~. - cfto*.. In the frame 
rotating at the drive frequency cj m , up to second order in the 
small parameters gj/Aj and gj/^j, this transformation yields 

2 ~ 

H eS = ArC^a ^^j-^-^j +Xj(T z j ^a 



Xj(jZ i ,„2-i2a> m t 



2 

8j* 



(a z e- lZ0Jmt + a n e i2aJmt ) + e m {a + a 1 ") 



(3) 



where A r = co r 



+Xj is the Lamb- shifted qubit 



frequency,^ = g 2 .{\/Aj + 1/2/) is the qubit state-dependent 
frequency shift and / = gigiiXl^j ~ 1/2/) is me interqubit 
coupling mediated by virtual excitation of the field. The terms 
that oscillate with exp(±i2oj m ) originate from the counter 
rotating terms acr~., a^cr* in the Rabi model. With the choice 
of co m « oj r for measurement on the qubits, we can safely 
ignore the qubit driving terms with amplitude gj€ m /Aj and 
gj€ m I^Lj lfT4l . In order to focus on the entanglement generated 
by measurement, we drop the term proportional to / since 
the possible measurement outcomes are the eigenstates of the 
interaction term cr\(T x v These measurement outcomes, on 
the contrary are not the eigenstates of the RWA's qubit-qubit 
interaction, which is proportional to cr+cr~ + h.c. lfl31 . 

In the Born-Markov approximation l25lL the dynamics of 
the system in the presence of dissipation and dephasing is 
described by a Lindblad form master equation j26l 

P = —IH^p] + KD[a]p + yijD[o-]]p + ?*L<D[o*]p 9 (4) 

where D[L]p = (2LpL^ - L^Lp - pL t L)/2. pit) is the density 
matrix of the total system, subjected to the loss of photons 
at rate k, energy relaxation of the qubit j at rate y\j and 
the dephasing of the qubit j at rate y^. Following l26l and 
neglecting the energy loss due to y\j, the intra-resonator field 
will evolve into a qubit- state dependent coherent states \a xy ), 
x, y = {g, e) with amplitude that satisfies 

a xy (t) = -ie m -i[A r a J y(t)+x*y(<x>y(t) + ^ 

(5) 

Herexxy = (xy\Xi°~\ +X20~ Z 2 \xy). These qubit-state dependent 
amplitudes will act as the pointers in the measurement where 
the information of the states of the qubits can be inferred from. 



In the situation where gj/co r «: 1 and Aj Zy, RWA is 
applicable to Eq. ([5]). In this regime, we can approximate 
Xj = g 2 j(l/ Aj+1 /I<j)by xJ WA = g^/A 7 and average out the fast 
oscillating exponential term of frequency 2oj m « 2oj r . This 
gives rise to the form 
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This result 



where ^ A = {mxr^i + *T "2' 
in the limit of RWA agrees with the expressions in Refs. 
lfT4l l26l . Comparing these two differential equations for 
the coherent states amplitude of the exact case [Eq. ([5])] 
and of the RWA case [Eq. (j6j>], we see that the former, 
apart from having an extra qubit state-dependent frequency 
shift of amplitude g 2 /^j, possesses an extra term involving 
the complex conjugate of the coherent state amplitude or* 
modulated by an exponential term with frequency 2oj m . This 
term originates from the two-photon process term in Eq. ([3]). 

To perform a parity measurement, the measurement must 
be able to distinguish the parity subspaces without destroying 
the quantum superpositions within the subspace. By choosing 
identical qubits (oj a j = a) a , j = 1,2) such that Aj = A, iZj = S, 
together with gj = g(xj = X) and A r = 0, we have*^ = ±2x 
for the states belonging to the even parity subspace *H+ and 
Xxy - for the states belonging to the odd parity subspace 
l27l . The corresponding phase space of the steady states 
is illustrated in Fig. [2] where the quadrature Q = lm[a xy (t)] is 
plotted against the in-phase component / = Re [0^(7)] of the 
transmitted field. 

In the homodyne detection of the resonator field, with the 
proper choice of local oscillator phase 0, one can measure the 
desired quadrature. The choice of = n/2, which corresponds 
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FIG. 2. (Color online) The phase space illustration of the steady 
state solution of a xy for both RWA (dashed lines) and exact (solid 
lines) with g/a) r = 0.001 (left) and g/a> r = 0.5 (right). The parity 
of the qubits can be inferred by choosing the relative phase of the 
local oscillator to be = tt/2, which project the field state to the Q 
quadrature. The parameters are A r = 0, e m = 0.5/c and g = 15/c with 
Q) r /27i = 2 GHz as fixed parameter. 
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to the Q quadrature, reveals only the information about the 
parity fT4l [T5l . However, there is still information in the 
quadrature orthogonal to the measurement, I (<p = 0). As 
shown in Fig. [2] the overlay between the coherent states 
amplitudes a ee and a gg in even subspace *H+ is not perfect. 
Hence this measurement is not optimal and leads to dephasing 
within the subspace <H+ . 

In Fig. [2] we compare the coherent states amplitudes 
in the regime of g/cj r = 0.001 and g/cj r = 0.5. Here 
we set ^//c = 15, the same value considered in a recent 
study of feasibility of error correction benchmarks for parity 
measurement l28l . To keep the validity of the dispersive 
regime, we choose g/A = 0.1 l29l and set e = 0.5/c to stay 
much below the critical photon number n cv a = (A/2g) 2 El. 
When g/oj r <?c 1, the contributions of the counter rotating 
terms are insignificant, hence the exact solution and the RWA 
one are indistinguishable [left column of Fig. |2). In the larger 
coupling strength regime where RWA is no longer applicable, 
the solutions differ from each other [right column of Fig. 
[2) . Deviations in both real and imaginary components of the 
coherent state amplitudes a ee are depicted in Fig. [3] In the 
exact case, the value of Re[a ee (t)] is closer to zero. Since 
Re[a ee (t)] = -Re[a gg (t)], this implies that the overlapping of 
the amplitudes increases. For the exact solution, there is also a 
slight increase in lm[a ee (t)], which means that the separation 
between parity subspaces is somehow larger in the exact 
solution. Meanwhile, due to the modulation of the oscillating 
term in the exact solution Eq. ([5]), small oscillations in the 
amplitudes persist even after the system reaches the steady 
state. 
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FIG. 3. (Color online) Dynamics of Re[a ee (t)] and lm[ar ee (0] vs time 
with coupling strength of g/a> r = 0.5. The other parameters are the 
same as in Fig. [2](Notice that the color code is different from Fig. [2]). 



III. ENTANGLEMENT GENERATION 

To facilitate the understanding of these distinctions between 
RWA solution and exact solution, we focus on one specific 
task, namely the generation of entanglement from the 
separable states with parity measurement. Consider the 



initial separable state (\g) + \e))l V2 <g> (\g) + \e))l V2 |0), 
which corresponds to the first qubit, second qubit and vacuum 
resonator field respectively. Upon the displacement of the 
driving field, the joint atom-field state can be written as 

m = 4= l<A+> K> + \\\ee) \a ee ) + \gg) \a gg )]. (7) 
V2 1 

The desired Bell states = Qeg) + \ge))/^2 and 

= (\ee) + \gg))l V2 can be obtained by Q quadrature 
homodyne measurement on the transmitted field. To quantify 
the distinguishability between the entangled Bell states, the 
average fidelity of the measurement is defined as 

F = P< H+ F (p+ ^P< H _F^, (8) 

where P<^ + and /Vl are me success detection probabilities 
for the even and odd subspaces respectively. The expressions 
for the fidelity of the Bell states and F^ + are given in 
Appendix A. In Fig. [4] we plot the average fidelity of the 
entangled states as a function of both g/K and g/a) r9 with 
a) r l2n fixed as 2 GHz. We observe that the average fidelity 




g/u r 



FIG. 4. (Color online) The average fidelity F of the parity 
measurement as a function of g/a> r and g/K with e = 0.5/c (top) and 
e - k (bottom) where oj r l2n = 2 GHz. The circle and square dots 
on the top figure (g/K =15) correspond to phase space in Fig. [2] 
The average fidelity of RWA (F RWA ) which corresponds to the value 
in the limit g/a> r — > is not shown since it does not capture any 
dependence in g/a> r and predicts a constant value. 
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improves as the ratio g/K increases. This is because with 
higher ratio of g/K, more information on the overall parity 
compare to the information on states within the subspaces is 
revealed fT4l . Notice that the maximum average fidelity for 
the choice e = 0.5/c is only about 0.84. A higher maximum 
average fidelity of about 0.97 can be achieved by choosing 
a larger driving field amplitude while remaining well below 
ftcrit, i.e. 6 = k. This is due to the fact that the coherent 
state amplitudes a eg and a ge of the odd parity subspace //_, 
which are proportional to the average number of photons (26), 
displace further away from the origin in the phase space as the 
photon number increases. Hence the separation between the 
parity subspaces increases, allowing the subspaces to be more 
distinguishable. 

Let us now follow the physics of the problem when g/a> r 
is varied. As the ratio g/co r further increases, as seen in the 
Fig. [4] the average fidelity F improves. This is due to the 
enhancement of the qubit state-dependent frequency shifty in 
the exact case, where the factor g/E = l/(A/g + 2aj r /g) in 
X increases as g/co r rises. This implies that the qubits shift 
the resonator frequency by a larger amount, thus decrease the 
interaction between the qubits and the driving field at bare 
resonator frequency oj r . Hence the dephasing in the even 
parity subspace decreases (HI . 

This enhancement in average fidelity is less significant for 
a higher ratio of g/K. For comparison we plot the average 
fidelity for g/K = 15 and g/K = 50 respectively in Fig. [5] 
(black and white dash line in the bottom column of Fig. [^J7 
The improvement in the average fidelity is significant for low 
ratio of g/K but not so much for a higher ratio of g/K. This 
is basically due to the fact that at large value of g/K, the 
overlapping between the states in the even subspaces is large 
and the parity measurement is close to optimal. As a result 
the advantage of the Rabi model over RWA in improving the 
average fidelity is less prominent, since the fidelity depends 
on the overlapping between the states in the subspaces and 
distinguisability between different parity subspaces. 



0.99 




0.2 . 0.3 

g/uj r 



FIG. 5. (Color online) Average fidelity for RWA solution (dash) and 
exact solution (solid) with g/K = 15 (red) and g/K = 50 (blue), which 
correspond to the white and black dash lines in Fig. [4] 



IV. CONCLUSION 

In summary, we have explored the possibility of the 
dispersive parity measurement in the ultrastrong coupling 
regime in circuit QED. We have shown that, in general, the 
fidelity of the parity measurement is enhanced in this regime, 
due to the additional frequency shift that depends on the state 
of the qubits. 
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Appendix: Average Fidelity 

To derive the average fidelity for the parity measurement 
on the initial state, we follow the procedure outlined in Ref. 
(30). Since we are measuring along the quadrature Q (cp = 
7r/2), the conditional state of the initial state Eq. ([7]) of the 
joint atom- field system for a measured p value may now be 
expressed as 



\V C (p))=^\* + )+^\ee)+^ 



V2 



, (A.l) 



where C xy (p) = G xy (p)K xy (p), with G xy = (2/;r) i/4 exp[-(/? - 
Im[of X} ;]) 2 ] and K xy = exp[-/Re[of X}; ](2/? - Im[of X}; ])]. After 
the homodyne detection of the resonator transmission, the 
joint system can be described by an unnormalized conditional 
density matrix p c (p) with the diagonal elements 

^ + \p c (p)\i(, + ) = \C eg (p)\ 2 /2, 

< < A_| P c (p)|^-> = o, 

(cf> ± \p c (p) |0 ± > = \C ee (p)\ 2 /4 ± l/8[C 2 ee (p) + Cf g (p)]. 

(A.2) 

To distinguish between the odd and even parity states we 
define the midpoint of the subspaces 



(Im[a^] + Im[a eg ])/2, 



(A.3) 



and assign the post measurement results to r H+ or <W_ if 
p > p m or p < p m . For the success probability can be 
explicitly shown to be 

P„ + = dpTr[p c (p)] 

= - (erfc ( <2{p m - lm[a eg \)) + erfc ( <2{p m - lm[a ee ]))) 



(A.4) 



similarly for P^, which is expected due to the symmetry 
between the subspaces over the quadrature Q. With this 
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success probability, the fidelity of obtaining Bell state |0+) 

from the *H+ subspace can be calculated as where b = Re[a ee ] and c = lm[a ee ]. For the orthogonal Bell 

state of from the subspace, the fidelity is found to be 

1 f°° 

F* + = — dp(cp + \p c (p)\cp + ) 

* < H+ Jp m y rPm 

1 (lm[a e J-lm[a ee ]\ ^ = \ dp (i// + \p c (p)\i// + ) 

= -erfc —- F( H J-oo 

4 \ V2 J 1 (Im[a eg ]-lm[a ee ]\ 

+ ie-^^M^^^-A) = rrfC i V5 1' (A ' 6) 
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erfc ( V2[ ^ m ^ gg ^ Im[ofg g ] ^ \\ ^ ^ With these expressions, the average fidelity F over the states 

\ \ - // |0+) and |0+) can be obtained by Eq. d8l). 
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